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I.  INTRODUCTION 


A.  GENERAL  DESCRIPTION  OF  THE  PROBLEM 


This  thesis  addressed  the  problem  of  estimating  the  distance 


which  must  be  traveled  to  visit  targets  which  are  randomly  distributed 


over  a large  area.  The  problem  is  interesting  because  many  practical 


search  operations  require  the  investigation  of  targets  after  detection. 


It  is  often  necessary  to  estimate  the  time  which  will  be  required  for 


the  entire  operation  before  the  targets  are  even  detected. 


Classical  search  theory  addresses  cases  where  it  can  be  assumed 


that  there  are  few  targets  in  the  area,  and  that  all  are  of  equal 


interest  when  detected.  Later  classification  may  rule  out  some  targets 


as  false  targets.  This  is  not  usually  obvious  to  the  searcher.  Any 


detection  will  stop  the  search.  Subsequent  classification  as  a false 


alarm  may  result  in  its  restart. 


A variation  on  this  case  occurs  when  there  are  many  false  targets. 


but  they  may  be  classified  immediately  upon  detection.  If  the  target 
may  be  immediately  classified  as  false,  then  the  search  can  proceed 


unimpeded.  False  targets  can  be  treated  as  if  they  had  never  been 


detected.  If  all  false  targets  are  immediately  classified  as  false. 


then  this  case  becomes  equivalent  to  the  classical  case  of  sparse 


targets.  The  search  continues  until  the  true  target  is  found  or  the 


area  is  exhausted.  When  the  first  target  is  detected  which  cannot  be 


dismissed  as  false,  the  search  is  stopped  and  that  target  is  pursued. 


There  is  a third  case  which  shall  be  addressed.  It  is  inter- 


mediate between  the  first  two.  It  can  be  defined  as  follows:  targets 


i 


can  be  easily  classified,  but  only  at  very  close  range.  The  searcher 
can  detect  targets  at  long  range,  but,  for  classification,  must  take 
a "close  look”  at  which  time  he  knows  Immediately  whether  or  not  he 
has  found  the  object  of  his  search.  If  he  has,  he  stops.  If  he  hasn't, 
he  keeps  searching.  If  the  targets  are  dense  enough,  or  the  searcher's 
detection  range  r^reat  enough,  the  searcher  may  always  have  uninvestigated 
targets  In  sight.  He  may  not  really  "search"  at  all  In  the  sense  of 
classical  theory;  he  may  Instead  travel  directly  from  target  to  target 
without  any  regular  "search  pattern."  The  problem  Is  now  one  of 
"Investigation  theory"  as  first  defined  by  Nuendorfer  Cll.  The 
Investigation  will  continue  until  the  desired  object  has  been  found, 
or  until  aTl  of  the  targets  have  been  visited.  In  either  event,  the 
Investigator's  track  will  have  gone  to  the  location  of  each  prior 
false  alarm. 

This  is  the  crux  of  the  problem.  If  the  positions  of  the  targets 
are  unknown  until  detection,  then  the  length  of  the  entire  path  will 
be  uncertain.  If  the  target  positions  are  random  variables,  then  the 
length  of  the  path  required  to  Investigate  the  targets  will  also  be 
a random  variable. 


B.  EXAMPLE 

The  case  described  above  Is  typical  of  many  Coast  Guard  operations. 
The  Coast  Guard  Is  often  required  to  search  for  a particular  boat. 
Usually  the  boat  can  be  positively  Identified  only  by  name  or  registra- 
tion number.  Small  boats  which  have  been  reported  overdue,  and  boats 
suspected  of  law  enforcement  violations  are  two  classes  which  must  be 
searched  for  In  this  manner.  There  may  be  many  boats  which  fit  a 
general  description,  and  the  searcher  may  have  to  check  out  all  of  them. 
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A simple  scenario  can  be  constructed  which  Illustrates  the  problem. 
Consider  a Coast  Guard  Cutter  with  embarked  helicopter  assigned  to 
patrol  a ten  thousand  square  mile  area  for  fishing  violations.  The 
helicopter  must  be  fueled  tonight  for  tomorrow's  first  light  patrol. 

The  aircraft  commander  must  decide  how  much  fuel  to  load.  Too  little 
fuel  means  an  Incomplete  patrol.  Too  much  fuel  means  that  something 
else  must  be  left  behind.  Experience  has  shown  that  rarely  are  more 
than  twelve  or  fewer  than  eight  ships  seen.  Ten  Is  the  average  number. 
The  aircraft  commander  has  no  way  of  predicting  the  future  position 
of  any  of  the  fishing  vessels,  but  he  must  still  make  an  Important 
operational  decision,  based  on  the  length  of  the  path  which  he  expects 
to  fly.  The  result  developed  in  this  paper  would  be  directly 
applicable  In  this  scenario. 
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II.  ANALYSIS 


A.  ANALYTICAL  FORMULATION  OF  THE  PROBLEM 

Two  basic  questions  must  be  addressed  first  off: 

1.  What  distribution  Is  to  be  assumed  for  the  positions  of  the 
targets? 

2.  In  what  order  are  the  targets  to  be  visited? 

In  this  thesis  It  Is  assumed  that  the  area  considered  Is  a portion 
of  an  infinite  Poisson  field  with  targets  distributed  with  spatial 
Poisson  density  p.  This  means  that  the  number  of  targets  In  any 
finite  area.  A,  will  be  a random  variable,  with  a Poisson  distribution 
and  mean  pA.  P.A.W.  Lewis  and  G.S.  Shedler  [2]  have  shown  that  the 
targets  occurring  in  an  area  as  the  result  of  this  process  will  be 
uniformly  distributed  within  the  area.  This  is  to  say  that  If  the 
entire  area  Is  divided  up  Into  many  small  boxes,  each  large  enough  to 
hold  only  one  target,  then  each  box  has  an  equal  chance  of  actually 
holding  one  target.  This  distributional  assumption  1s  attractive 
because  It  assumes  the  minimum  amount  of  Information  about  the  target's 
location.  In  this  sense  it  Is  the  "worst  case";  more  Information  can 
only  Improve  the  situation.  In  actual  Coast  Guard  practice,  there  are 
many  searches  conducted  when  uniformity  Is  the  best  assumption  which 
can  be  made;  the  actual  distribution  of  targets,  including  false  targets, 
cannot  be  distinguished  from  what  would  have  been  expected  from  a 
uniform  distribution. 

A related  question  Is  that  of  target  motion.  In  this  thesis  It 
will  be  assumed  that  the  targets  move  so  much  more  slowly  than  the 
Investigator  that  they  may  be  considered  stationary. 
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The  question  of  the  order  In  which  targets  are  to  be  visited  is 
more  complicated.  If  ♦•he  targets  are  stationary,  there  is  one  unique 
optimum  order  resulting  in  the  minimum  total  path  length.  The  problem 
of  determining  this  optimum  ordering  is  called  the  "travelling 
salesman  problem."  The  travelling  salesman  solution,  while  ideal, 
is  not  practical  for  large  problems.  There  is  no  efficient  method  of 
solution.  Every  current  method  is  fundamentally  equivalent  to  trying 
all  possible  orders  [3].  There  are  some  shortcuts  for  small  numbers 
of  points;  the  fact  that  points  on  the  perimeter  will  be  taken  in  the 
relative  order  that  they  are  found  around  the  perimeter  is  an  example. 
But  these  shortcuts  do  not  work  for  large  numbers  of  points  with 
complex  geometry.  Determining  the  optimum  path  through  twenty  points 
will  require  a tremendous'  amount  of  calculation. 

This  thesis  employs  the  myopic  decision  rule  in  which  the  investi- 
gator will  always  go  to  the  nearest  point  which  has  not  yet  been 
visited.  This  decision  rule  is  justified  on  two  grounds: 

1.  It  is  easily  employed  in  actual  searches;  the  searcher  need 
only  remember  what  targets  have  been  investigated,  and  then  select  the 
closest  of  the  remainder. 

2.  Balut  [4]  has  shown  simulation  results  indicating  that  the 
path  resulting  form  the  myopic  decision  rule  is  only  about  152  longer 
than  the  optimal  path  for  investigations  where  the  investigator  is 
much  faster  than  the  targets.  This  procedure  is  often  used  in  actual 
searches. 

Finally,  the  assumption  is  made  that  the  investigator's  detection 
range  is  much  greater  than  the  inter- target  spacing.  This  means  that 
the  investigator  will  always  know  the  location  of  the  nearest  uninvesti- 
gated target. 
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B.  MODEL  FOR  AN  UNBOUNDED  MYOPIC  INVESTIGATION 

Because  the  area  considered  is  a portion  of  an  infinite  Poisson 
field,  it  is  profitable  to  first  consider  the  properties  of  an  investi- 
gation conducted  in  an  infinite  field.  The  investigator  will  start 
at  one  of  the  points  in  the  field  and  then  visit  its  closest  neighbor, 
then  the  closest  neighbor  remaining,  and  so  on.  The  investigator  will 
always  travel  to  the  nearest  uninvestigated  target  (n.u.t.).  Real 
searches  never  detect  an  infinite  number  of  targets,  but  this  is  not 
a bad  approximation  to  the  case  when  the  area  is  very  large  and  there 
are  many  more  targets  detected  than  can  be  investigated  in  the  course 
of  a mission. 

Consider  a Poisson  field  with  spatial  density  p.  Call  the  set  of 
points  I = {Pg,  P-j.  Pg.  P3.  •••.  Pp>  an  investigation  starting  at  Pg.  If 
p^  = Pg,  i.e.  the  investigator  returns  to  the  point  where  he  started, 
then  the  investigation  is  a tour.  The  investigation  is  conducted 
according  to  the  following  rule:  the  investigator  will  always  visit 
the  closest  point  to  which  he  has  not  yet  been.  Let  r^  be  the  distance 

k n 

from  p^_^  to  p^.  Let  r^..  Let  t^  = r^.,  when  p„  = Pg. 

is  the  cumulative  distance  convered  in  visitina  the  first  k ooints; 

t is  the  length  of  a tour  visiting  all  the  points  and  returning  to 
n 

its  starting  point. 

The  expected  length  of  the  first  radius,  r^,  may  be  calculated  as 
follows.  The  probability  of  a target  being  found  in  a small  area,  dA,  is 

P(dA)  = pdA. 

Thus  the  probability  that  a target  will  not  be  found  is 

f(dA)  » 1 - pdA. 
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Let  r(A)  be  the  probability  that  there  Is  no  point  In  some  large  area  A 
The  probability  that  there  Is  no  point  In  (A  + dA)  Is 

T{^)  • RdA)  - r(A)(l  - pdA). 

But  this  Is  also  T'(A  dA) 

r(A)(l  - pdA)  - r(A  + dA) 


where  c Is  the  arbitrary  constant  of  Integration 


The  area  considered  so  far  has  been  perfectly  general,  of  arbitrary 


shape. 

Now  consider  A to  be  the  area  circumscribed  by  a circle  of  radius  r 
2 

A * irr  . Let  A^  be  the  smallest  such  circle  tangent  to  an  adjacent 
point.  Let  r-j  be  its  radius.  The  tail  distribution  of  A^ , (t)  is 

defined  as  follows: 


where  it  Is  the  area  circumscribed  by  a circle  of  radius  x 


Define 


The  expected  value  of  r, , E(r. ),  is  defined  as 


Integrating  by  parts  yields 


where  S’  (x)  » P(r,  > x) 

’^1  ' 

p » Spatial  Poisson  density 


Thus  the  expected  radius  of  the  smallest  circle  containing  an  adjacent 
point  Is 


E(r,) 


1 

2*^ 


Figure  II-l 

The  expected  values  of  successive  jumps  are  much  more  difficult 
to  determine.  The  distribution  In  the  fresh  (unsearched)  area  Is  still 

T(A)  » 

but  the  new  area  is  no  longer  a circle.  Let  x^  be  the  actual  distance 

between  Pq  and  Its  nearest  neighbor,  p-j.  Then  a circle  of  radius  x-j, 

centered  at  Pq,  Is  known  to  contain  no  other  points  besides  Pg.  The 

new  point,  p.j,  lies  on  to  this  circle.  Area  Included  In  this  circle 

"doesn't  count"  In  the  search  for  new  points  and  must  be  Ignored  In 

calculating  the  distribution  function  for  the  next  Interpoint  distance, 

5!.  (x).  Let  H.(x)  be  the  function  whose  value  ir  the  new  (uninvestigated) 
1 

area  lying  within  a circle  of  radius  x,  centered  at  P^. 
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Thus 


j f 

i t 


Hn(x)  « irx^  and  (x)  * e 

□ 


•pHn(x) 


It  can  easily  be  seen  that 


H-jCx)  ■ it(x^  - x^^)  if  (x  >_2x^) 


^ /T 


i(x  - X,  ) - (2x,(^)cos  - 2J^ 


4x 


If  (x  < 2xi)  where  x^  Is  the  actual  observed  distance  between  Pq  and  p-j 


Let  (x)  » PCr,  > x). 
r2  £ 


Then,  as  before 


M • e 

'^2 


•pH^(x) 


and 


E(r2) 


-pH,(u) 
e du 


/, 


2x, 


-1  X 


exp{-p{ir(x  - x^  ) - [2x^( — - l)cos  ^ ^ 


2x, 


/: 


4x, 


]}}dx 


r " 2 2 

I exp  { -pir(x  - x^  )}dx 
J2x^ 


This  expression  Is  difficult  to  evaluate  analytically,  and  the  expected 


values  of  r^,  r^  ...  are  even  more  difficult  to  obtain.  This  Is  because 


the  expected  value  for  rj^  depends  not  only  upon  r^  and  rg  ...  » 

but  also  upon  the  angles  between  successive  radii.  An  Investigation 
which  doubled  back  on  Itself,  Figure  II-3(a),  would  have  a different  next 
expected  radius  than  one  which  did  not. 
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Figure  I I -3 


While  the  total  sanitized  area  is  the  same  in  both  cases,  it  is 
not  distributed  according  to  a common  function  of  radius. 

Because  the  distributions  and  moments  are  so  difficult,  approxima- 
tions would  be  useful.  An  argument  can  be  advanced  to  support  the 
conjecture  that  E(r^),  i = 1,  2,  3,  ...  are  asymptotically  linear  as 
i -»■  ».  Because  new  points  will  be  located  in  an  un investigated  area, 
the  process  will  tend  to  creep  away  from  sanitized  areas  and  toward 
unsanitized  areas.  This  will  usually  put  the  investigator  at  the 
boundary  of  a fresh  area,  with  a constantly  expanding  sanitized  area 
behind  him.  Asymptotic  linearity  may  be  shown  easily  in  a restricted 
model  based  on  this  behavior.  Consider  a completely  deterministic 
process  where  points  appear  on  a straight  line.  The  spacing  between 
the  points  will  be  determined  as  follows: 
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1..  win  be  arbitrary, 

2.  r^  Will  be  such  that 

irr^^  » “ 


^ ( 


11 

I 

' ! 
j 

ii 


2r 


1}cos 


-1  *^1 


2r 


1-1 


1-1 


2r 


1-1 


/: 


4»*i-i 


for  all  In  1 > 1 . 

This  Is  equivalent  to  circumscribing  a circle  of  radius  r^  about  the 

first  point,  and  then  adding  other  points  so  that  the  new  area  added 

to  the  process  at  each  step  Is  equal  to  the  Initial  area  enclosed  In 

the  first  circle.  The  process  can  be  visualized  as  a sequence  of 

circles  proceeding  down  a line.  Each  new  circle  Is  centered  on  the 

Intersection  of  Its  predecessor  and  the  centerline.  Each  circle  Is 

sized  so  that  Its  area,  less  the  area  of  Its  predecessors.  Is  equal  to 
2 

r^  . Refer  to  Figure  I I -4. 


Numerical  calculation  for  all  Initial  values  selected  between 
r-j  » 1 X 10“^  and  r^  ■ 1.0  x 10^®  have  shown  that  the  sequence 
*"1*  *^3*  "^n*  •**  asymptotic  to  1.2814  r^. 


I 
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While  one  Mould  never  expect  to  see  an  actual  realization  develop 
along  a perfectly  straight  line,  an  extensive  computer  simulation  by 
Laldlaw  [5]  does  show  that  many  Investigations  do  develop  In  a way  that 
could  be  characterized  as  random  deviations  on  either  side  of  a mean 
line  of  progress.  It  seems  Reasonable  that.  If  this  were  so,  and  the 
Investigation  did  proceed  as  variations  about  a mean  line,  then  the 
actual  limiting  value  might  be  r^  cr^,  where  c ^1.2814.  Laldlaw's 
simulation  yielded  c-j  » 1.405  based  on  a myopic  investigation  conducted 
in  a finite  plane  containing  10,000  points,  with  the  investigation 
stopping  when  the  boundary  of  the  plane  was  closer  than  the  nearest 
target,  and  with  the  sweep  width  much  greater  than  the  expected  target 
spacing. 

The  assumptions  for  this  result  were: 

1.  The  points  visited  are  located  in  an  infinite  Poisson  field. 

2.  Detection  range  is  much  greater  than  target  spacing. 

3.  Myopic  search  is  employed. 

Under  these  assumptions  the  length  of  the  total  path  appears  to 
be  asymptotically  convergent  to 

* ncr, 
n 1 

, n(1.405) 

2«^ 

The  model  is  asymptotically  linear  in  n and  is  called  the  asymptotic 
linear  model  of  the  unbounded  investigation. 


19 


C.  APPLICATION  TO  THE  BOUNDED  INVESTIGATION 


I 

t I 

\ i 


I 

} 

L 


The  linear  asymptotic  model  cannot  be  directly  applied  to  the 
bounded  Investigation.  The  linear  model  requires  that  the  spatial 
Poisson  paramenter,  p,  remain  constant  In  the  unsearched  area.  This 
trill  not  be  true  In  a finite  area,  because  the  finite  area  will 
contain  only  a finite  number  of  points. 

Nevertheless,  the  linear  asymptotic  model  does  suggest  an  Indirect 
approach.  Consider  a bounded  Investigation  to  consist  of  two  parts: 

1.  An  Investigation  of  the  points  near  the  center  of  the  area 
where  the  unbounded  linear  model  may  be  applied. 

2.  An  Investigation  of  the  remaining  points  near  the  boundary. 

It  seems  reasonable  to  suppose  that  the  bounded  Investigation  might 
first  behave  like  the  unbounded  case,  when  there  are  a large  number 
of  points  from  which  to  draw.  After  a while  the  points  will  have 
been  "thinned  out"  considerably,  and  the  Investigator's  path  will  take 
him  back  to  "mop  up"  the  remainder. 

A simulation  was  constructed  to  test  this  conjecture.  A number 
of  points  were  generated  in  a square.  An  Investigation  was  started 
at  a point  near  the  center.  The  path  took  the  Investigator  to  all 
the  points  In  the  area,  and  then  back  to  his  Initial  point.  Such  a 
closed  path  Is  called  a tour.  Plots  of  t|^/kE(r^)  - vs  - k for 
k * 1,  2,  3 ...  n,  n * 10,  30,  60  shown  In  Figure  III-5.  A cumulative 
plot  of  t^/nE(ri)  - vs  - n Is  shown  In  Figure  III-6. 

Two  things  can  be  seen  In  these  graphs: 

1.  The  "plateau"  In  the  center  of  the  curves  In  Figure  III-5  becomes 
more  pronounced  as  more  points  are  added  to  the  area.  This  Is  a 
reasonable  result. 
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KEY:  o*#=20,  -•-^=40,  »-Ht=60 
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Figure  II-5 
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2.  The  final  cumulative  length  is  approximately 

t^  - 1.95^-!!|^  » 1.95nE(r^).  n > 7.‘ 

Plots  of  repeated  simulation  trials  suggest  that  t^  is  normally 
distributed  with  mean 

T 1 nc 

\ * T 

c 2 n nc  ''A 

S.  » 0.06  ^ 

Stemleaf  plots  for  the  cases  n * 20,  40,  and  60  are  shown  in 
Figure  III-7. 

This  is  the  expression  for  t^  when  n is  known.  Generally  n will 
not  be  known,  but  the  assumption  that  the  search  area  is  taken  from 
a Poisson  field  will  enable  the  calculation  of  E(t^)  from  p and  A. 
The  expected  number  of  points  occurring  in  an  area  A is  E(n)  * pA. 
Thus 

E{t^|n)  » 1.95uE(r^) 

1 


* 1.95  n 


2/? 


and 


E(t„)  - i 1.95  k-!— Mk) 
" k-0  2/T  " 


where  f„(k)  ■ P(ri  • k) 


« 1.95  — ^ 

w 

I 1 

• 2/T 

k«0 

• 1.95  — J- 

E(n) 

2/T 

* 1 .95  — 2- 

• pA 

2/~p 
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similarly, 

« 0.06  £=5— 

III.  EXAMPLE 

The  procedure  developed  may  be  applied  to  the  example  given  In 
Section  II. 

Assume  that  the  search  area  Is  a square  area,  100  miles  on  a side, 
with  a spatial  Poisson  parameter  of 

p * 1.1  X 10"^. 

The  expected  length  of  a tour  will  be  given  by  * 1.95 

thus 

E(t^)  * 323  miles. 

To  check  this  result,  eleven  points  were  selected  at  random,  with 


one  being  the  searcher. 

The  positions. 

measured  from  the  southwest 

comer,  were: 

Index 

£ 

L 

1 

32.1 

51.1 

2 

60.5 

00.4 

3 

90.5 

17.3 

4 

49.8 

18.2 

5 

18.4 

99.2 

6 

65.3 

72.9 

7 

40.8 

12.8 

8 

51.6 

22.2 

9 

69.7 

99.3 

10 

58.0 

74.0 

0 

18.0 

19.8.  (searcher) 
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These  points  are  plotted  In  Figure  III-l.  It  Is  assumed  that  the 
searcher's  detection  range  Is  greater  than  the  greatest  Inter- target 
spacing. 

The  myopic  Investigation  starting  at  point  "0"  Is  shown  by  the 
solid  line  In  Figure  III-2.  I » {0.  7,  4.  8,  2,  3,  6,  10,  9.  5.  1.  0}. 
The  actual  length  of  the  myopic  Investigation  starting  at  pt  0 Is 
328.66.  In  this  case  the  myopic  path  Is  also  the  optimal  path  through 
all  the  points,  and  328.66  Is  the  minimum  length.  In  this  case  n » 11 
and 

S.  ^ « .06  ^ 1650 


S.  « 40.62 
^n 

So  the  actual  value  Is  well  within  one  standard  deviation  of  the 
predicted  value. 
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I Point  Locations  for  Example  Problem 
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IV.  SUMMARY  AND  CONCLUSIONS 


This  thesis  has  presented  a method  of  estimating  the  length  of  a 
path  required  to  visit  all  of  the  targets  In  a search  area  under  three 
conditions: 

1.  The  targets  are  distributed  uniformly  within  the  area. 

2.  The  "nearest  uninvestigated  target"  (myopic)  decision  rule  is 
used. 

3.  Detection  radius  Is  much  greater  than  Inter-target  spacing. 
Simulations  results  have  shown  that  the  expected  distance  of  a complete 
tour  of  all  points  under  these  assumptions  is  predicted  very  well  by 
a simple  function  of  the  area  and  number  of  targets.  If  that  number 
Is  known,  or  by  the  area  and  average  target  density.  If  the  number  of 
targets  Is  not  known. 
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